Thermoelectric power in one-dimensional Hubbard model 
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The thermoelectric power 5* is studied within the one-dimensional Hubbard model using the linear response 
theory and the numerical exact-diagonalization method for small systems. While both the diagonal and off- 
diagonal dynamical correlation functions of particle and energy current are singular within the model even at 
temperature T > 0, S behaves regularly as a function of frequency lo and T. Dependence on the electron 
density n below the half-filling reveals a change of sign of 5 at no = 0.73 ± 0.07 due to strong correlations, in 
the whole T range considered. Approaching half-filling S is hole-like and can become large for U ^ t although 
decreasing with T. 

PACS numbers: 71.27.+a, 71.10.Fd, 72.15.Jf 



INTRODUCTION 

The thermoelectric power (TEP) or the Seebeck coefficient 
S can show quite anomalous behavior in metals with strongly 
correlated electrons. This has been recognized very clearly 
in the normal-state of high-Tc cuprates, where numerous ex- 
perimental studies 1 1, 2] reveal several unusual features of the 
TEP S{T): a) In the underdoped regime with small concen- 
tration of holes n?j <C 1 introduced into the reference antifer- 
romagnet (AFM) S is positive and large, showing a decrease 
with temperature dS/dT < 0, except at very low T ^ Tc 
where the usual Fermi-liquid (FL) behavior S oc T sets in, b) 
as a function of hole doping the TEP S also decreases and ap- 
proaches zero or even changes sign in the overdoped regime. 
The latter phenomena can be attributed to the effect of strong 
correlations, since for noninteracting (or weakly interacting) 
electrons one would expect negative and small S < Q within 
the same regime. 

Although one of the easiest transport quantities to measure, 
the TEP proved to be a hard challenge for a theoretical con- 
sideration and understanding. Apparently, it is even more 
demanding than some other transport properties as, e.g., the 
electrical conductivity cr, which are also anomalous in many 
correlated systems, and in cuprates in particular. Since the 
Boltzmann transport theory can fail completely for strongly 
correlated electrons, the general linear response approach has 
to be followed for 513]. Nontrivial expressions have been ob- 
tained in this way in the limit of high temperature |4] where 
S{T cx)) approaches a finite value. The limiting value 
of iS* is also known for the extended Hubbard model with an 
infinite-repulsion i7 — > oo f?] and for more more compli- 
cated models |5]. An extension of S{T) to finite T < oo 
has been studied within the U = oo Hubbard model using the 
retraceable-path approximation f3]. A numerical analysis of 
S{T) within the planar t-J model has obtained several fea- 
tures relevant to cuprates |7]. The latter calculation is based 
on an observation that the ratio of off-diagonal and diagonal 
dynamical correlation functions between the energy and par- 
ticle current, respectively, is nearly constant. Using the exact 
diagonalization of small systems the TEP has been studied 



also for models with the orbital degeneracy j|8|]. Another line 
of approach is to express S{T) in terms of spectral functions. 
This is well posed approximation within the dynamical mean- 
field theory |9, 10], valid for the Hubbard model in the limit 
of large dimensions. Using spectral functions obtained via the 
fluctuation-exchange approximation the latter approach has 
been applied also to the 2D model 1 1 1 ] and further improved 
by the inclusion of vertex corrections tl2ll . 

The one-dimensional (ID) Hubbard model, which is the 
subject of this study, is known to possess special features rela- 
tive to systems with higher _D > 1. As ID models of interact- 
ing electrons in general, it complies to the phenomenology of 
Luttinger liquids and exhibits the spin-charge separation I isil . 
Using the Bethe-Ansatz solution for low-energy spinon and 
holon excitation branches, the low-T behavior of S has been 
considered in Ref.| 14]. Results appear quite unusual since S 
is claimed to remain positive at all electron densities n < 1 
below half filling. This study indicates that even S{T 0) 
is an open question within the ID Hubbard model, and even 
more the behavior at higher T. On the other hand, it has been 
recognized that the integrability of the ID Hubbard model im- 
plies singular transport quantities [15, 16, 17]. In particular, 
it has been shown that away from half-filling, i.e., for n 1, 
the charge and energy (heat) transport remain dissipationless 
even at T > 0, i.e., in the d.c. limit the electrical as well as 
the thermal conductivity both diverge, o-q ^ oo and Aq ^ oo 
flTil . Since S{T) is the ratio of two correlation functions, it is 
expected to remain well defined and finite 1 15]. 

In order to clarify the behavior of the TEP within the ID 
Hubbard model, we present in the following results of the nu- 
merical study of the dynamical iS'(u;, T), based on the linear 
response theory, generalized to allow for singular transport 
quantities. Dynamical response functions are evaluated nu- 
merically for finite chains with L sites using the approach of 
exact diagonalization. In particular, for the largest systems 
we apply the finite-temperature Lanczos method (FTLM) 
which enables us to consider a model with up to L = 14 sites . 
We concentrate on the regime of strong correlations U > 2t. 

The paper is organized as follows: In Sec. 2 we present 
the general linear response formalism for transport quantities 
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within the ID Hubbard model, with the special care for singu- 
lar behavior in the d.c. limit lu ^ 0. Different limits are also 
considered. Results of the numerical study of finite systems 
are presented in Sec. 3 and conclusions are given in Sec. 4. 



LINEAR RESPONSE THEORY 



It follows from Eq.(|6j that - 



-Ht/L and 



TnE = -j- + H.C.) 
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Xl{(c/+i,sC/s + clci-i,s)ni-, + H.C.}, (8) 



We consider the ID Hubbard model 
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(1) 



with the nearest neighbour (n.n.) hopping term Ht and the 
onsite repulsion [/ on a chain of L sites. The particle current 
has the usual form (we set further h — 1 and the lattice spacing 

ao = 1), 



Jn = t ^(«c|^i_sQs + H.C.), 

l,s 



(2) 



whereas the ener gy c urrent can be expressed in terms of n.n. 
energy operators ilTlflsIl 



Je = i^[hi,hi+i], H ^y^ hi. 
I I 

It then follows O from Eq. Q 

Je - -i^ ^(ic;Vi,sC/-i,. + H.C.) + 



(3) 
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-T + 4^ci-i^s)nL-s + H.C.}. (4) 
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Within the linear response formalism we calculate general 
transport coefficients Lij |3], defined by Ji = ^jLijFj. 
Here, i, j = n, Q, where Jq ^ Je — fJ-Jn is the heat current, 
/i is the chemical potential, and related generalized forces are 
given by F„ = -{Vfi + eo£)/T, Fq = -VT/T^. Since the 
transport coefficients in the ID Hubbard model are singular 
in the d.c. limit lITsillTll . it is essential to generalize the for- 
malism to finite frequencies, dealing with response functions 
Lii ( (jj). If one introduces appropriate polarization operators 
fl?,"2tfl Pr, = Ini, Pe - E; and Pq = Pe - y.Pn 
(valid for open boundary conditions), we can write within the 
linear response theory 
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where lo^ 



and 



x..M = -T / '^^e- *([J.(i),J,(0)]) 



(5) 
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while r„Q = TnE — t-'-Tnn- Note that the dynamical electri- 
cal conductivity, more frequently discussed in this context, is 
given by a{uj) ^ e^i„„(tj)/T. 

Finally, we get for the dynamical TEP jsl 



S{u;) 



1 Lnqjuj) ^ 
eoT L„„(tj) 

1 f {tue) - XnE{^) 
gqT V (r„„) 



(9) 



In most model systems (and in real materials) we are dealing 
with the situation where transport coefficients at T > are 
not singular within the d.c. limit, i.e., Lij{uj 0) — L^^ are 
finite. This is related to the ergodic behavior of the system and 
Eq.(|3 in this case requires the equality of the d.c. adiabatic 
and isothermal susceptibilities, Xiji^ 0) = xjj = {Tij) 
d. 

It has been shown that in this respect the ID Hubbard model 
is singular Jvh . as well as some other ID integrable quantum 
many -body models f\3\. Let us define the generalized stiff- 
ness 
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Within the ID Hubbard model -Djj(T > 0) remain finite 
away from half-filling n ^ 1 |17]. The general argument 
has been given with the finite overlap of current operators 
{JiQn) where Qn,n — are conserved quantities due 
to integrability. Note that in a macroscopic system there is 
an infinite number of operators Qn- In particular, Q3 corre- 
sponds closely to Je, obtained by the replacement in Eq.0 
U/2 U. In spite of singular Lij{u! 0) the d.c. TEP 
S = S{llj 0) is expected to remain well defined and finite. 
From the representation chosen in Eq.([^ we can express S 
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Dr. 



eoT \D„ 



(11) 



There are some limits which give a nontrivial test to our re- 
sults presented below. For noninteracting electrons at T > 
it is straightforward to evaluate S. Since J„ and Je are con- 
served quantities (constants of motion), we get from Eq.(0 
that Xnn(w) = XnE{uj) = and 



Tnn = j fie)g{e)ede, 
TnE = J f{e)gie) {2e^ - 4f) de, (12) 
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where .g(e) = 2{4t^ -e^)-^^"^ /n, f = l/(exp[(e-/i)/T] + l) 
and fj, is fixed by the density n ^ J f{e)g{e)de. 

In the limit of large UbT » {U, t) the result coincides with 
that of noninteracting electrons at ksT ^ t, 

5 = _^ln^. (13) 

eo n 

On the other hand, a nontrivial generalized Heikes formula 
is obtained for the case of large {ksT^ U) ^ t but UbT <C U 

i, 

5=-^ln^ii^. (14) 

eo n 

It has been recently observed |7] that in the most frus- 
trated regime of strongly correlated systems the transport co- 
efficients LnE and L„„ are closely related. Numerical results 
for the planar t-J model thus reveal (at the intermediate dop- 
ing of the reference antiferromagnet) that in a broad range of 
Lo,T the ratio R{ijo) — i„£;(ti>)/L„„((jj) is nearly constant. 
Extending this behavior to T ^ a simple relation for S{T) 
would follow, which is very attractive for application in com- 
plicated correlated systems |21|. It is of interest to which ex- 
tent within the ID Hubbai'd model R{T) ^ Aifi/^Drm be- 
haves similarly. 

RESULTS 

We analyse dynamical transport coefficients Lnn{^) and 
Lnsi'-^) and consequently S{lo) by performing the exact di- 
agonalization of the model Eq. ([0 on a finite chain with L 
sites and with periodic boundary conditions. Taking into ac- 
count the number of electrons and the total spin projection 
as conserved quantities as well as the translation symme- 
try (wavevector q) systems up to L = 10 sites (all fillings) 
can be fully diagonalized to get results at arbitrary T. Using 
the FTLM one can analyse larger systems at T > jQ and 
in the following we present results obtained for L — 12, 14 
with the number of basis states within a symmetry sector up 
to Nst ^ 800.000. All Ne < L (for > L we employ 
the electron-hole symmetry of the model) are investigated in 
order to perform the grand-canonical averaging which allows 
for a continuous variation of n — N^/L. In the present FTLM 
application we use typically M = 100 Lanczos steps and ran- 
dom averaging over R ^ 20 samples within each symmetry 
sector. 

Within FTLM, as well as in any exact calculation of a small 
system, it is important to realize that results have a restricted 
thermodynamic validity due to finite-size effects, which begin 
to dominate at low T < Tfg- Clearly, Tfs depends on the 
size of the system, but as well on the density and character of 
low-T excitations. It is characteristic that ksTfa/t ^ 0.2 (at 
n ^ 1) m the ID Hubbard chain with L = 14 is larger that 
in a corresponding 2D (square lattice) Hubbard model with 
N = 4x4 = lOsites |22] (where it was found fcsT/s/t - 0.1 



at the intermediate doping of the reference insulator) in spite 
of an evidently shorter linear extension. The difference can be 
explained with the anomalous behavior of the planar Hubbard 
model (in the strong-correlation regime U ^ t)\22] and of re- 
lated 2D i- J model 1 7] which show a non-FL behavior close to 
half-filling, i.e., a large degeneracy of low-lying states results 
in a large entropy and small Tfs- On the other hand, the ID 
Hubbard model at low T follows the phenomenology of Lut- 
tinger liquids ITill . with well defined excitation branches of 
spinons and holons, resulting in an increased but not anoma- 
lous entropy. 




FIG. 1: Transport coefficients Re L„„(lj), Re I/„_B(tj) and TEP 
[^(lj)!, calculated within the ID Hubbard model with U/t = 8, 
L = 14, T = 0.5 t, for different electron densities n close to half- 
filling. Results obtained via FTLM are broadened with S = 0.25 t. 

Let us first discuss the frequency-dependent L„„, LnE and 
S. In Fig. 1 we present Re L„„(ti>) oc cr. Re LnE{'-^) as well 
as |S'(ti')|, calculated using the FTLM in a i = 14 system 
for the typical strong-correlation regime U = 8 t for fixed 
ksT = 0.5 t and different n close to half-filling. Spectra are 
smoothed with an additional broade ning 5 = 0.25 t. As al- 
ready found in previous studies ClIlJ, finite D„niT > 0) 
(away from half-filling n ^ 1) is reflected in a singular (delta- 
function) contribution to Re £„„(aj — > 0) in Fig. 1. The 
incoherent part pronounced for uj > At, mainly due to the 
transition from the lower to upper Hubbard band, decreases 



away from n ~ 1, consistent with the free electron propaga- 
tion when approaching the empty band, n ^ 0. An analogous 
structure is observed in the off-diagonal Lnsi^), whereby the 
singular contribution DnE can be negative as well, and in fact 
shows the change of sign upon decreasing n. The resulting 
(absolute value) 15(^^)1 evaluated from Eq.(|9} is presented in 
Fig. 1. In the strong-correlation regime it seems to be char- 
acteristic (close to half-filling) that S{u!) shows a pronounced 
resonance at uj ^ ujq with ujq linked to the pseudo-gap in the 
incoherent part. The resonance is determined by the mini- 
mum of \Lnn{i^o)\, i.e., by the vanishing Im i„„(wo) = 
whereby also Re L„„((jJo) is small. This results in a sharp 
peak in |S'(w)|, with a simultaneous change of phase in com- 
plex S{uj). The same effect we find for U = At. Away from 
w ~ the dynamical TEP |5'(a;)| behaves quite smoothly 
approaching both (different) nontrivial limits u; = Q and 



■ U/t=2 




FIG. 2: TEP S (in units fcs/eo) vs. T for noninteracting electrons 
for different densities n, calculated from the analytical result <12lllt . 
compared to the result for a small system L — 14. 

In Fig. 2 we present the reference d.c. TEP S{T) for nonin- 
teracting fermions, evaluated from Eq. J12I1 11 1 for different n. 
We also display the result of a grand-canonical calculation of 
S" on a finite-size system with L = \A sites. The deviations of 
the latter from exact results at low T can serve as an estimate 
for Tfs, relevant for FTLM calculations presented below. It 
should be noted that in general at T ^ small-system results 
for S do not get a meaningful limit, in particular they do not 
reproduce correct FL S* oc T behavior It follows from Fig. 2 
that ksTfs ^ 0.2 — 0.3 t, depending on n. Note also that 
5 = for n = 1 due to the particle-hole symmetry of the 
model, which remains valid also for general [/ ^ Ii23i. It 
should also be observed from Fig. 2 that the FL regime with 
5 oc T is restricted to quite low ksT < 0.4 t even for nonin- 
teracting fermions, in spite of much broader band W = At. 

FTLM results for S{T) for different U/t = 2, 4, 8 and n = 
0.3 — 0.9 are presented in Fig. 3. For a given system of L sites 
and fixed T we calculate using FTLM S{lij 0) values for all 
sectors 1 < < L (taking into account the symmetry of the 
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FIG. 3: S vs. T for the ID Hubbard model with U jt = 2,4,8, 
respectively, as calculated for different n with the FTLM for L — 
12, 14 chains. 



model for L+1 < < 2L — 1) and perform grand-canonical 
averaging to get a continuous variation S{T, n). We show in 
Fig. 3 results for two systems, L — 12, 14, which allow to 
estimate Tfs for C/ > cases. 

Several observations can be made on the basis of S{T) re- 
sults in Fig. 3: a) For n < 0.6 we get for all considered U/t 
S{T) < 0, with values essentially identical to that of nonin- 
teracting electrons in Fig. 2. b) The effect of [/ > becomes 
well pronounced close to half-filling, where the TEP becomes 
hole-like, i.e. S > 0, and at the same time large. At large 
U ^ t, the TEP S{T) stays constant in a large T window, 
whereas for smaller U, e.g. U/t = 2, it is suppressed already 
for modest ksT ^ t to the noninteracting value S ^ 0. c) It 
is rather difficult to reach the FL regime where S oc T. Still, 
away from half-filling, i.e. for n < 0.6, numerical results ap- 
pear consistent with the behavior S{T 0) 0. This is, 
however, not evident for n ~ 1, in particular not for largest 
U ^8t. 

In Fig. 4 we display the alternative representation of the 
TEP results S{n) (for L — lA system only) for various 
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n 



FIG. 4: S vs. n for U/t = 2, 4, 8 and different T, as calculated with 
FTLM for L = 14 chains. For U/t = S case we present also the 
Heikes value, Eq. <14> . 



ksT/t = 0.3,0.5,1.0 and same U/t as in Fig. 3. Again, 
for n < 0.6 the TEP is electron-like S < and behaves 
as for noninteracting electrons. The interesting part is for 
n > 0.6 where U > induces a transition to a hole-like 
5 > for n > uq. In general, no is dependent both on T 
and U. The limiting value, reached for U ^ fc^T ^ i, is 
given by the generalized Heikes formula, Ea.( ll4> . with the 
value no — 2/3. Moreover, as seen in Fig. 4, uq remains 
nearly constant, i.e. no ^ 0.7, in a broad range of T in the 
strong-correlation regime U/t > 4. Only for a weak repul- 
sion U/t — 2 case the crossing starts to approach the nonin- 
teracting no = 1. Since at n = 1 we get S = for all T, 
there must be a maximum S* = S{n = n*) within the hole- 
like regime. The location n* as well as the value S* depends 
strongly on U and T. The limit, Ea. (ll4> . gives n* = 1 with 
divergent S* — > oo. For U — 8t our results for ksT < 0.5 t 
show a good overall agreement with the latter limiting behav- 
ior Note that here due to large charge gap at U = 8 t, the 
maximum appears very close to half-filling, i.e., n* ^ 1 with 
S* ^ ks/eo. With decreasing U/t as well as with increasing 
T the maximum becomes less pronounced and moves away 
from n* = 1. 



CONCLUSIONS 

Above we have presented results for the dynamical and d.c. 
TEP within the ID Hubbard model, as calculated by using the 
FTLM for small systems with up to L = 14 sites. Several 
comments and conclusions are in order: 

a) The main restriction in the validity of results comes from 
the use of small systems. Finite-size effects in S and other 
quantities start to dominate results for T < Tfs- Although 
the FTLM algorithm involves also the random sampling over 
initial configurations fT], it does not essentially increase the 
numerical error of S and only slightly increases Tfs- 

b) It is quite characteristic that observed Tfs is weakly depen- 
dent on [/ and n, i.e., ksTfs/t - 0.2-0.3forL = 14. This is 
strikingly different from the Hubbard model (or t-J model) in 
2D, where Tfs varies substantially with concentration n in the 
strong-correlation regime |7, 22], i.e., Tfs has a minimum at 
the 'optimum' doping where the quantum antiferromagnet is 
frustrated by mobile holes. The latter regime is characterized 
by the large entropy, hence also low Tfs- In the ID Hubbard 
model there seems to be no such phenomenon, consistent with 
the low density of excitations (spinons and holons) within the 
Luttinger liquid. 

c) We have investigated only the regime below half-filling 
n < 1, since due to the particle-hole symmetry of the model, 
Eq. (□, it follows S{2 - n) = -S{n)- 

d) Independent of U at low filling n < 0.6 the TEP is electron- 
like, S < 0, and essentially equal to that of noninteracting 
electrons in the whole regime of T considered. Although the 
result appears plausible, it contradicts the conclusion based on 
the exact Bethe-Ansatz solutions 1 14] that for large U > t the 
TEP is hole-like S" > for all n < 1. 

e) TEP becomes strongly dependent on U as well as on T 
when approaching the half-filling. In particular, it changes the 
sign at no, which in a broad regime of U, T stays close to the 
Heikes result no ~ 2/3, but increases towards no — > 1 with 
the vanishing ?7 0. S{n) reaches a maximum at n = n*, 
located in the regime no < n* < 1. Both n* and S* are 
dependent on U and T. In particular, S* can become very 
large in the strong-correlation regime U ^ t where n* 1 
and S{n < n*) follows Eq. ( flU down to low T. 

f) The low-T variation of the TEP is expected to follow the 
standard FL behavior S oc T- It appears quite difficult to 
enter and investigate the latter regime using small-system di- 
agonalization. The linear dependence is restricted to a rather 
narrow window ksT < 0.5 t even for free fermions, as seen 
in Fig. 2. For U ^ t close to the half-filling this window is 
expected to become even narrower, since the relevant energy 
scale is that of spinons with the characteristic exchange cou- 
pling J = At'^/U- So it is not surprising that our results in 
Fig. 3 for largest J7 = 8tatn~ldo not reveal the onset of 
the linear regime, down to ksTfs ~ 0.3 t- In any case, the 
low-T behavior of the TEP remains an open question which 
could be possibly settled by an analytical calculation fT?]. 

g) It is a relevant question to what extent the behavior within 
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ID Hubbard model reflects the general features of transport 
and of TEP in higher dimensional systems with strongly cor- 
related electrons. Of particular interest is the relation with the 
physics of quasi-2D high-T^cuprates and with their anoma- 
lous electronic properties llj, yj, which have been the mo- 
tivation of the intensive research of microscopic models of 
correlated electrons. Apart from results based on the approx- 
imation of infinite dimensions |9, 10], there are few model 
calculations of S in 2D systems |7] and the question of the 
variation of iS'(r) in prototype models is not settled. Still, 
some differences and similarities between ID and higher D 
systems seem to be evident. Transport in ID Hubbard model 
is specific due to its integrability, showing in the divergence 
of d.c. coefficients (w 0). Nevertheless, the TEP S is 
well behaved and even some deviation from integrability (e.g., 
by adding additional next n.n. hopping in H) should not in- 
fluence its behavior dramatically. The difference between ID 
and 2D (or higher D) models is in the nature of low-lying ex- 
citations. While in the ID model the latter follow the scenario 
of the Luttinger liquid, the situation in 2D is more involved 
and not fully understood yet. At least numerical investiga- 
tion reveal a very degenerate quantum state close to optimum 
doping ItI E^I . One consequence appears to be nearly lo and 
T-independent ratio R = LnE/Lnn 01 • In contrast, in ID 
model R = DnE /Dnn remains rather T-dependent (roughly 
as much as /i), so S cannot be well expressed solely in terms 
of/x(T). 

h) Nevertheless, there are important similarities between 
the TEP in ID Hubbard model and available 2D results [7J 
and moreover experimental results within the normal state 
of cuprates Ql |^ . In the strong correlation regime C/ 3> t, 
S becomes hole-like for n < 1. In a broad range of T 
the change of sign appears at no — 0.73 ± 0.07, close to 
the limiting Heikes value n = 2/3. At the same time the 
maximum value S* > can become large and remains 
as such down to quite low T, giving the TEP in strongly 
correlated system also the potential applicability. At the same 
time, in the hole-like regime the T dependence of the TEP is 
quite pronounced. In particular, dS/dT < at T > T*, also 
consistent with experiments on cuprates, although T* appears 



higher than in 2D systems 0] ■ 
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